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ABSTRACT 



We analyze the velocity gradient across the jets of the young stellar object 
(YSO) NGC 1333 IRAS 4A2 and explain it by decomposing the two opposite 
jets to two opposite sub-jets. The two sub-jets have the same radial velocity 
and the angle between them is only several degrees. Each of the two sub-jets is 
composed of two components. We also show that the alternative interpretation 
of jets' rotation is unlikely to account for the velocity gradient. The two sub- 
jets directions are constant, and there is no indication for precession. The line 
connecting the centers of the sub-jets defines a preferred direction in a plane 
parallel to the accretion disk launching the jets. We suggest that the preferred 
constant direction in the accretion disk is determined by the semi-major axis of a 
highly eccentric orbit of a brown dwarf or a massive planet companion. At each 
periastron passage the companion perturbs the accretion disk with a quadrupole 
mode, e.g., two spiral arms, that lead to the launching of the two sub-jets. We 
predict that a careful monitoring of the system will reveal a periodic activity 
with a period of few months to several years. 



INTRODUCTION 



Asymmetry in the line of sight velocity across some jets launched by young stellar objects 
(YSO s) has been inte r preted as caused by a la r ge scale rotation of the material in the jets 
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2011a ). The se measurements, howeve r, are quite c hallen ging and compli- 
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cated (e. g.. ICoffev et al.ll2012l). In two e arlier pap ers ([Soker 



2005 



results of Bacciotti I fl2002[ ). ICoffev et al.l (120041 ) and boffev et al.l ( 120071 ) and argued that the 
observations do not support, and even contradict, the interpretation of jet s rotating around 



their symmetry axes. The claims for jets' rotation in the YSOs HH 211 ( ILee et al.l ( 120071 ). 
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and Ori-S6 (IZapata et al.l |2010| ) were found to be quite problematic as well (ISoker 



Fast rotation around the jet axis is predicted by the magneto-centrifugal ac celeration 
(MCA) model for jet launching (e.g., lAnderson et al.l l2003t iFerreira et al.l 120061 ) . In this 
model the magnetic fields that are anchored into the accretion disk-star system play a dom- 
inate role in accelerating the jet's material from the accretion disk by a lever-arm torque. 
In this model the specific angular momentum carried by the jet's material is ~ (ARf) 2 Qf, 
where Rf is the radius at the foot-point of the jet on the disk, Qf is the angular velocity at 
Rf, and A > 1 is the lever- arm. 



Instead of the rotation interpretation, ISoker I (120051 ) proposed that interaction of the jets 
with a twisted-tilted (wrapped) accretion disk can form the observed asymmetry in the jets' 
line of sight velocity profiles. The jets interact with the ambient gas residing on the two sides 
of the disk, e.g., a weak disk outflow or corona. Such a model can explain the asymmetry in 
the line of sight velocity across the jets, avoiding the need to invoke jet rotation with a huge 
amount of angular momentum. 

Another mo del that can account fo r such observations without invoking jets' rotation 
was proposed by ICerqueira et al.l (120061 ) who assumed a precessing jet whose ejection ve- 
locity changes periodically with a period equals to the precession period. Practically, the 
dependance of the jet's expansi on velocity o n direction around the symmetry axis leads to 



the same effect as the model of ISoker I (120051 ). 



In a recent paper IChoi et al.l (j2011al ) argue that the velocity gradient across the jets of 
the YSO NGC 133 3 IRAS 4A2 can be interpreted as jets' rotation, and that this supports the 
disk- wind models (jPudritz et al.ll2007l ). In section [2] we summarize the properties of the out- 
flow from NGC 1333 IRAS 4A2, and in section [3] we present our two sub-jets interpretation. 
Our short summary is in section HI 



2. THE OUTFLOW FROM NGC 1333 IRAS 4A2 



The following properties of the jets of NGC 1333 IRAS 4A2 (hereafter IRAS 4A2) are 



(Choi et al. 


2010). 


Choi et al. 


(2011a) 



of M acc = 1.6 ± 0.9 x 1O- 6 M yi 
velocit y of the jet t hey o bserve in SiO is as that of the parts observed in molecular hydrogen 
line by lChoi et al.l (120061 ): 71 km s -1 in the plane of the sky. From that and the assumption 
that the velocity gradient is due to rotation they derive the specific angular momentum of 
the gas in the jets to be lj ~ 1 — 2 x 10 21 cm 2 s _1 . Assuming the magnetohydrodynamic 
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wind model for the rotating jets (lAnderson et al.ll2003l ; iPudritz et al.l 120071 ) they calculate 
the jets' foot-point to be at Rf ~ 2AU. At the location of this foot-point the specific angular 
momentum of the disk material is 



l d (R f ) = 1.8 x 10 19 
and its specific kinetic energy is 



0.08M M 



1/2 



R 



f 



2 AU 



1/2 



cm 2 s 1 ~ 



0.01 lj, 



e d = 1.8 x 10 



ii 



O.O8M 



Ri 



2 AU 



erg ~ 3.6 x 10 e, 



■o-> 



(2) 



where ej ~ (70 km s ) /2 is the specific kinetic energy of the gas in the two jets. 



The above relations have two implications. First, the 'lever-arm' is A = [Ij/l^Rf)] 1 ^ 2 a 
10 to A = [tj/ediRf)] 1 / 2 ~ 17. The lever -arm is defined as the ratio of the Alfven radius 
inward to which the outflowing gas is controlled by the disk's magnetic field, to the foot-point 



radiu s. Although values of A > 10 were suggested in the literature (e..g, ICasse &: Ferreira 



2000l . who by 'lever-arm' refer to a different quantity) and cannot be ruled out, this is a 



very large value for A which has a typical value of 2 < A < 3 (e.g., lAnderson et al.l 12003 
Pudritz et al.l 120071 ). We also note that the outflow velocity of Vj ~ 70 km s -1 is about the 
escape velocity from a YSO of mass M* = O.O8M , suggesting to us a jet origin very close 
to the star. 

Second, angular momentum and energy conservation with the above ratios imply that 
the jets mass outflow rate is limited to Mj < 0.0036M acc ~ 6 x 10~ 9 M Q yr _1 . There is 



no rel iable derivation of the mass loss from IRAS 4A2. S caling the results of iBlake et al. 
( 119951 ) to a distance of 235 pc a nd an age of ~ 50, 00 yr (jChoi et al.l 1201 lal ). gives a mass 



outflow rate of ~ 1O~ 6 M yr~ x . iLefloch et al.l (119981 ) finds the number density in the wings 
of high- velocity components of the outflow, v ~ 30 km s -1 , to be ~ 0.6 — 2 x 10 6 cm -3 close 
to the star. For jets' radius of 1" the calculated mass loss rate is ~ 10 _5 M o yr -1 . This 
seems to be too high. The conclusion here is that although better determination of the mass 
outflow rate is required, available data, with their large uncertainties, do suggest that the 
mass outflow rate is larger than the mass outflow rate allowed in the rotating jet model. 



Soker I (120051 ) proposed that interaction of the jets with a twisted-tilted (wrapped) ac- 



cretion disk can form the observed asymmetry in the jets line of sight velocity profiles. This 
type of interaction predicts that the western edge of the northern jet be symmetric to the 
eastern edge of the southern jet. Indeed, as seen in Figure [Uthe eastern edge of the northern 
jet (left side of upper jet in the figure) and the western edge of the southern jet (right side of 
the lower jet in the figure) have wiggling boundaries, while the western edge of t he northern 
jet and the eastern edge of the southern jet have relatively sharp boundaries (jchoi 2005 
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Choi et al.l 12011 ai) . This not only shows the type of symmetry expected from the model of 



Soker I (120051 ). but also that an interaction with the surrounding gas might account for the 



velocity gradient perpendicular to the jets' axis. In the next section we examine in greater 
detail the jets, and propose yet another possible process that operates in the launching 
process of the IRAS 4A2 jets. 



A DOUBLE-JET MODEL FOR THE JETS OF NGC 1333 IRAS 4A2 



Choi et al.l (j2011al ) used the Very Large Array to observe the bipolar jets of the protostar 
IRAS 4A2 in the SiO t> = 0J=l— »0 line. They investigated the kinematics of the two 
opposite SiO jets by placing slits perpendicular to the jets' axis; these were terme d cuts. The 



j ets an d the positions of the cuts (marked by numbers) are shown in Figure [TJ IChoi et al. 



fcOllaf ) produced position- velocity (PV) diagrams for each cut, wh ere they pos itioned the 
zero displacement of each cut at the emission peak according to IChoil ( 120051 ). The PV 
diag ram for each of th e perpendicular slits is presented in Figure [2] here, taken from Figure 
2 of IChoi et al.l (1201 laf ). where the colored crosses are our addition which we discuss below. 



Noticing that the sense of velocity gradient is the same in all cuts they argued for 
jets' rotation. To demonstrate their interpretation that the velocity gradient is due to jets' 
rotation, they drew a velocity gradient for each cut; these are the black dotted lines on the 
PV diagrams of Figures |2j We criticized this interpretation in section [2j We hence aim at a 
different interpretation of the velocity structure. Examining the PV diagram of the different 
cuts, we notice that there is a multi-peak structure in every one of them. (1) There are two 
clear separated structures in cuts 10 and 11. (2) There are three or four peaks in most cuts. 
(3) In cut 8 there are only three peaks, but the central one is the strongest among all peaks. 
This is as if two peaks overlap each other. 

We try to fit the peaks with the following structure. We assume the existence of two 
sub-jets, each composted of two components. Here a sub-jet refers to both northern and 
southern parts. Namely, a sub-jet has two opposite parts relative to the central star. The 
same holds for a 'component'. We model each of the four components as a straight line of gas 
on the two opposite sides of the star, as if it is a jet of its own. These four components are 
drawn on Figure [H The axis of each component is parametrized by two angles. First is the 
inclination relative to the plane of the sky, 9. For 9 > the northern part points away from 
us, i.e., red-shifted. Second is the angle on the p lane of the sky rela tive to component la. The 
values of 9 and <j) are given in Tabled] Based on lChoi et al.l ( 120061 ) we take all components to 
have the same radial expansion velocity of Vj = 71 km s _1 . We compute the PV diagram of 
the four components and draw them as colored crosses on Figure [2] where the velocity and 
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position is at the center of the cross. The crosses themselves represent our crude estimations 
of the errors in the derived values of 9 and <fi. From our trials of different fittings we crudely 
estimate the errors to be of few degrees. 

The four components do not perfectly fit the observed PV diagram as can be seen in 
Figure [21 First, there is a substantial interaction of the jets with the ambient gas. This is 



clear ly evident from the sharp bend of the northern jet which occurs further out (IChoi et al. 



20061 ; not seen here). Other regions also show some interaction with the surroundings, as 
evident from the wiggling boundaries of the jets (see discussion in section [2] here). Second, 
this interaction or jittering of the jets' source might cause random change in the expansion 
directions of the different components. 

We notice that we can substantially improve the fitting by displacing components la and 
lb together in cut 8, and components lie and lid in cut 9. The displacement is A8 = —3.7°, 
namely redward in the southern jet, for components la and lb in cut 8, and A6 = 2.0°, 
namely blueward in the southern jet, for components Ic and Id in cut 9. The PV diagram 
with these displacements is presented in Figure [3J The displacement in cut 8 explains why 
the peak there is as twice as high as in the other cuts. As the simultaneous displacement of 
components la and lb accounts quite well for the PV diagram in cut 8, we identify them as 
a single sub-jet I. As well, we identify components Ilc+IId as sub-jet II. 

Our main conclusions from the simple fitting, although not perfect, are as follows. 



1. The two opposite jets of IRAS 4A2 are actually composed of two bipolar sub-jets. Each 
of the bipolar sub-jets is composed of two components. 

2. We can make the fitting with two sub-jets that have the same radial velocity, and with 
a very small angle between them. These suggest a common launching site. 

3. There is no sign of precession. The sub-jets and the two components composing each 
sub-jet maintain an average constant direction. 



Table 1: The angles of each component 



Sub-jet 


Component 


9 





I 


la 


12.8° 


0° 


I 


lb 


11.2° 


5.0 


II 


lie 


4.6° 


6.5 


II 


lid 


7.1° 


7.7 
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4. The structure of four-components and two sub-jets exists on both sides of the star 
(north and south jets). This shows that the sub-jets structure does not result from a 
random collision and stochastic interaction with the surrounding gas. 

5. The relatively large displacement of one sub-jet we identify in cuts 8 and 9, according 
to our interpretation, can result from either change in launching direction or from 
interaction with the surrounding gas. The interaction with the surrounding gas can 
explain the small displacements seen in the different cuts and in small changes in 
structure and intensities. However, the large displacement found in cuts 8 and 9 are 
more likely to result from instabilities and stochastic behavior of the launching process 
that lead to a temporary change in the velocity (or launching angle) of one sub-jet. As 
well, the existence of the two sub-jets also points to a complicated launching process. 
We speculate on a possible cause in the next section. 

6. No extra component of jet rotation is required in the two sub-jets interpretation. 

We note that the sub-jets structure is found in the PV diagram and not in the image. 
This is possible because the sub-jets are almost in the plane of the sky, i.e., sin 9 <C 1. For 
higher values of sin# > 0.5 the contribution of the radial velocity to the Doppler shift is 
higher, and changes in the radial velocities of the sub-jets as they expand might overwhelm 
the differences between the sub-jets and components. 

4. DISCUSSION AND SUMMARY 

Our interpretation of two sub-jets structure as presented in the previous section requires 
that the launching process has a preferred direction that does not change with time, at least 
not for the time period span by the jets of ~ 500 yr. This excludes the possibility that the 
preferred direction is the direction toward a companion in a circular orbit. This is because 
for a slowly varying direction toward the companion the orbital period must be thousands of 
years or more. Such a companion will be at a too large distance to influence the launching 
process of the jets. A close companion changes its direction on a short time scale. 

However, if the orbit is eccentric the semi-major axis of the orbit defines a preferred 
direction that can maintain a constant direction for a very long time. We propose, therefore, 
that the launching process of the jets in IRAS 4A2 is substantially influenced by the presence 
of a companion on a highly eccentric orbit. 



Ardila et al.l ( 120051 ) present the interaction of a disk with a close binary star on a 
parabolic orbit. This passage creates two large-scale spiral arms in the disk. The spiral 
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arms extend inward to a distance of ~ 0.1r p , where r p is the periastron distance, and pre- 
serve a more or less constant direction for most of the time. The dense-spiral arms structure 
in the gas component of the disk lasts for a time of ~ 0.5P p , where P p is the orbital period 
of a circular orbit of radius r p . The two spiral arms can lead to the formation of the two 
sub-jets we identify in the outflow from IRAS 4A2. 

We can estimate plausible parameters for the orbit. Taking the IRAS 4A2 jets' velocity 
of 70 km s _1 to be the escape speed fro m their launching radius r^, we find tl — QRq, where 



the mass of the star is M* = O.O8M ( jChoi et al.ll2011bl ). Let the companion periastron be 
then at r p ~ 3 — 10r^ ~ 20 — 60i? Q , and the apastron be at £ ^> 1 times this distance. The 
orbital time of the binary system (assuming the companion to have a mass of M2 << Mi = 
O.O8M ) is 

J./T \ -1/2 / \ 3/2 /-, . ( \ 3/2 

M* \ f r p \ ( 1 + 4 x 



Each jets' launching episode lasts for a few months. Over the hundreds of years life time of 
the jets, the many launching episodes will be smeared to continuous jets. While the system 
is near apastron, the disk near the primary star can rebuild itself. 

One prediction of the model we proposed is that a careful monitoring of the system 
will reveal a periodic activity with a period time between few months and several years. 
Observations should be "lucky" enough to catch the system near a periastron passage, when 
massive launching takes place. The companion is a brown dwarf or a massive planet. 

We cannot elaborate on how the two spiral arms in the disk lead to the formation of two 
sub-jets. Simply, there is no agreed-upon model for the launching of jets that we can use to 
derive the process. As for the two components of each sub-jet, numerical simulations of the 
gravitational interaction of the companion with the disk around the primary star should be 
carried out, probably in three-dimensions. Such simulations can reveal the presence of higher 
harmonics in the accretion disk in addition to the spiral arms, as well as the development of 
instabilities. But again, the way the structure in the disk influences the jets formation is a 
more complicated process. 

The main claim of this paper is that the velocity gradient across the jets of the YSO 
IRAS 4A2 cannot be accounted for by jets' rotation (section [2]). They are much better 
interpreted as arising from a sub-jets structure (section [3]). A similar interpretation can 
hold for the outfl ow HH 797 located in the IC 348 cluster in Perseus. In a recent paper 



Pech et al.l (120121 ) interpret the velocity gradient across the jet of this outflow as evidence 
for jet rotation. We suggest that a careful examination might reveal a sub-jets structure in 
that outflow. However, in the case of HH 797 interaction with the ambient gas seems to play 
a major role, and interpretation is more complicated. In any case, a two sub-jets structure 
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is clearly seen in their PV diagrams. The origin of the two sub-jets structure might be the 
presence of a close companion on a highly eccentric orbit. 

This research was supported by the Asher Fund for Space Research and the E. and J. 
Bishop Research Fund at the Technion. 
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Fig. 1. — The color composi te image of the NG C 1333 IRAS 4A2 bipolar jet in the SiO 
v = J = 1 — > line from IChoi et al.l (j2011al ) . The white lines acros s the je ts that are 
labelled by numbers, the 'cuts', are from the original figure of IChoi et al.l (j2011al ). while the 
four lines along the two opposite jets are our addition. These four lines mark the two opposite 
sub-jets (J and //), each composed of two components that we identify in the sub-jets (see 
text). 



Position Offset 



Fig. 2. — The position- velocity (PV) diagrams from lChoi et al.l (j2011al ). The colored crosses 
are our addition. The center of each cross represents the Doppler shift and position of a sub- 
jet component along the appropriate cut. Each of the four symbols of the crosses represents 
one component, as marked on two of the panels. 
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Fig. 3. — Like Figure El but with a displacement of sub-jet I in cut 8 and of sub-jet II in cut 
9, as marked by the vertical black arrows (see text). 



